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The spectral emission of neutrinos from a plasma in nuclear statistical equilibrium (NSE) is inves- 
tigated. Particular attention is paid to the possible emission of high energy (>10 MeV) neutrinos 
or antineutrinos. A newly developed numerical approach for describing the abundances of nuclei 
in NSE is presented. Neutrino emission spectra, resulting from general Fuller, Fowler, Newman 
(FFN) conditions, are analyzed. Regions of T-p-Ye space favoring detectability are selected. The 
importance of critical Ye values with zero net rate of neutronization (Ye) is discussed. Results are 
provided for the processing of matter under conditions typical for thermonuclear and core-collapse 
supernovae, pre-supernova stars, and neutron star mergers. 
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I. INTRODUCTION 



Neutrino cooling is ofparamount importance in the 
modern astrophysics [H, 0, 0] . It governs late stages of 
stellar evolution, especially massive stars (>8-10 M©) 
[1, [1] , red giant cores B, white dwarfs [3] , core-collapse 
supernovae [sL lol. ITol. [Tiril^. ITsj and (proto)neutron stars 
[14| . Neutrino emission is irnportant in mergers involving 
neutron star [l^, [l^ [13, [H, , the dense accretion disks 
of Gamma Ray Bursts (GRB) modelsH HH, [H, type 
la supernovae 23] and X-ray flashes [24l |. 

Usually, neutrinos carry away energy, and only the to- 
tal neutrino emissivity, i.e. amount of energy carried out 
by neutrinos is of interest. The neutronization induced 
by the net j/g — flux is crucial for understanding of 
the nucleosynthesis. Therefore previous research on NSE 
neutrino emission [1^ focused on: (1) Ug — Vf. particle 
emission rates and (2) total ~\- energy carried out 
by the neutrinos. We would like to extend this analysis 
to cover spectral/flavor properties of the NSE neutrino 
flux. 

In known research a detailed treatment of the neutrino 
emission is done for core-collapse simulations [1^ [l^l . On 
the other hand, it is frequently neglected for other astro- 
physical objects (e.g. la supernovae). Nowadays more 
interest is dedicated towards spectral properties of the 
neutrino flux. The neutrino energy is important for core- 
collapse supernovae, for the neutrino-induced nucleosyn- 
thesis (i^-process, [11, [2^, [s^l ) , neutrino oscillations, and 
for the detection of neutrinos in terrestrial experiments. 
The last area is poorly explored. The neutrino spectrum 
for neutrino cooling processes rarely is treated in rigorous 
way. Typical procedure is to use some more or less jus- 
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tified analytic forms for the neutrino energy spectrum. 
There are parameters that are found from known neu- 
trino emissivity and the average neutrino energy. In this 
paper we continue our former investigation [3ll . [32,] to 
flnd spectral properties for important neutrino emission 
processes, We proceed now to processes involving weak 
nuclear (5 transitions. 

Neutrino cooling processes can be separated into two 
classes. There are (1) thermal processes including e~e^ 
pair annihilation, massive in-medium photon & plasmon 
decay and neutrino photoproduction, and (2) weak nu- 
clear processes (i.e. decays and captures). We 
would like to point out that for all thermal processes 
(pair, plasma, photoproduction, bremsstrahlung, neu- 
trino de-excitation of the nuclei) the neutronization rate 
vanishes, i.e. the change of the proton/neutron ratio is 
due exclusively to weak nuclear processes. Class (1) pro- 
duces all flavors while (2) only v,. and i>e. However, neu- 
trino oscillations can mix flavors. Information on thermal 
and weak components might be destroyed. It happens 
somewhere between emission and interaction/detection. 

We assume that matter is transparent to neutrinos. 
Therefore, weak nuclear processes often tend to domi- 
nate neutrino emission of hot and very dense plasma. In 
particular, electron captures by both protons and heavy 
nuclei are progressively more intense. With growing den- 
sity, the Fermi energy Ep ~ can become larger than 
the capture threshold (Q-value) , for increasing number of 
nuclear species. High temperature additionally enhances 
emission. Many of the nuclei remain in the thermally ex- 
cited states. Matrix elements for these weak transitions 
are often large. For temperatures above ~0.5 MeV, a sig- 
nificant fraction of equilibrium positrons builds up. This 
causes a strong flux due to e"*" captures, particularly 
on free neutrons. 

In contrast to thermal processes, determined entirely 
(including energy spectrum) by the local thermodynamic 
properties of matter (e.g. temperature kT and electron 
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chemical potential He), weak nuclear processes depend 
also on abundances of nuclei. This renders the task of 
calculating neutrino spectrum difficult to achieve. This 
is especially true for evolutionary advanced objectS"'^. All 
that we can say for rapidly evolving object, is that the 
neutrino spectrum emitted from plasma is of the form: 



P 



nipAk 



(1) 



Here ipk represent (assumed known, from theory or exper- 
iment) spectral shape of single nuclei neutrino emission, 
and Xk{t) set of usually unknown and rapidly varying 
abundances. Tracking of the required required number 
of a few hundred abundances is possible at most in sim- 
plest one-dimensional models, to our knowledge. 

Fortunately, if the temperature becomes high enough, 
nuclei begin to ,,melt" due to photo-disintegrations. Nu- 
clei re-arrange due to strong interactions into the most 
probable state favored by the thermodynamics [s^]. This 
is the Nuclear Statistical Equilibrium (thereafter NSE) 
approximation [sl]. The timescale req uired to achieve 
NSE is temperature-dependent [H, [s^. It can be ap- 
proximated as [st! ]: 



0.2^179. 7/Tc,-40. 5 
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(2) 



where p is the density in g/cm^, Tg ~ T/IO^K and T 
is the temperature in K. Eq. ^ provides one of the 
most important constraints limiting the use of the NSE 
approach. We assume implicitly in ^ that Ye ~ 0.5 
[35, _36] . Therefore in a plasma with the value of Ye which 
is far from 0.5 caution is required. Both under- an over- 
estimate is possible. The timescale is of the order of the 
age of the universe, tnse ~ 10^ years, for kT = 0.2 MeV 
and Tnse ~ 10~^ seconds for kT = 1 MeV. In the core 
of a typical pre-supernova star with p = 10^ g/cm"^ and 
kT = 0.32 MeV we have tnse — 2 days. A typical du- 
ration of the Si burning stages depends on stellar mass 
and varies from few hours to 3 weeks. During the ther- 
monuclear explosion of type la supernova in the flame 
region temperatures grow up to kT — 0.4 . . . 0.6 MeV, 
the timescale tnse — 5 milliseconds, and the explosion 
time is of the order of 1 second. 

The weak transmutation rate between protons and 
neutrons is denoted by Ye, 



Y 



dYejt) 
dt 
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where: 



Ak 



If the hydrodynamic timescale is longer than tnse and 
Ye change slowly^ then we can safely assume a qua- 
sistatic evolution in the three-parameter space. Usu- 
ally'^ these parameters are temperature kT, the den- 
sity p and the electron fraction Ye- For the given triad 
{kT, p. Ye) we are able to determine abundances of all nu- 
clei. This approximation is widely used in " iron" cores of 
pre-supernova stars, supernovae, nuclear networks, ther- 
monuclear flames and nucleosynthesis studies. Under 
NSE conditions the neutrino emission is not much differ- 
ent from thermal processes (especially if Ye — 0), and no 
prior knowledge of abundances is required. This allows 
e.g. for post-processing of models with a known history 
of temperature, density and electron fraction. If Ye{t) is 
not known we still can use the NSE approximation as- 
suming some value, e.g. Ye = 0.5 for symmetric nuclear 
matter. The composition (and therefore neutrino emis- 
sion) is extremely sensitive to small changes in Ye in the 
most interesting range of Ye = 0.35 . . . 0.55 and relatively 
low temperatures of kT < 0.5 MeV. One method to over- 
come this problem is to use the so-called tracer particles 
built into simulation to remember the thermodynamic 
history of matter. In the next step one then finds the 
history of Ye- Another application of the NSE neutrino 
emission, described in [23, is the subgrid-scale model of 
nuclear flame energetics in thermonuclear supernovae. 

This paper is organized as follows: in Sect. [TT] we 
discuss spectra of individual nuclei under conditions of 
high temperature and density. We use solar ^Be neu- 
trinos as an example. In Sect. IIIII we use NSE and 
get neutrino emissivities and energy spectra, using FFN 
[381 [sol. |40|. [4]| weak rates. Final section comprises con- 
cluding remarks and a programme for future theoretical 
neutrino astronomy (calculations of the neutrino spectra 
and so on). 

For details related to the implementation of NSE the 
reader is directed to the accompanying paper, that is 
submitted to Atomic Data and Nuclear Data Tables [421 . 



II. NEUTRINO SPECTRUM FROM /3 
PROCESSES IN THERMAL BATH 

Bahcall [l^, laid fundamental theoretical founda- 
tions in the context of Solar neutrino spectrum. Later 
work is upgrade for the results of Bahcall concerning the 
number of nuclei involved, better nuclear data etc. With 
notable exception of the Sun and geo-neutrinos [i^ 
a rigorous treatment of the neutrino spectra from indi- 
vidual nuclei is usually ignored in astroparticle physics. 
Core-collapse simulations use parameterized approach, 



^ This situation is however very difficult to describe using statisti- 
cal methods. Variety of astrophysical objects and processes make 
it closer to complex systems rather than gases. 



^ Slow in the sense of eq. not actual weak rates A^, Xa^ , which 

may be very hig h. 
^ As noted by 1251 . relativistically invariant triad T—ns — Ye where 

riB is conserved baryon number density may be used if General 

Relativity formulation is required. 
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cf. e.g. [23, |43|. Unfortunately, in the case of multi- 
peaked neutrino spectrum this approach simply does not 
work, cf. Fig. 1 and related comments in [43]. The an- 
tineutrino spectrum is computed only for free neutrons, 
in applications known to the author. 

The spectrum of neutrinos emitted from single nuclei 
in astrophysical plasma depends strongly on the tem- 
perature and the chemical potential of electrons (and 
positrons if kT ^ nie = 0.511 MeV or larger). The tem- 
perature is large in typical evolutionary advanced astro- 
physical objects (pre-supernova or supernova, for exam- 
ple). We will study neutrino spectrum in this regime. On 
the contrary, for the solar interior, kT — 1.35 x 10~^, /ie — 
0, and this makes little change with respect to laboratory 
experiments. 

Let us begin with typical example of the continuum 
electron capture process: 



Normalized y,, spectrum 
0.5 



^Bc 



Li 



We make assumptions concerning the infinite nucleus 
mass and we neglect various correction factors (screen- 
ing. Coulomb factor). Then the capture rate is pro- 
portional to the constant matrix element multiplied by 
the so-called phase space factor $: 



£l[£,-^Q)^{£,-/^QY-m 
l-hexp[(£^- AQ-^)/fcr] 



-e(£:,-AQ-me), 
(4) 

where E^, denotes the neutrino energy for e capture and 
£1, is the antineutrino energy for capture. AQ is the 
energy difference between initial and final states (both 
can be excited) and m^. is the electron rest mass. The 
chemical potential /ie of the electron includes rue, and 
therefore for positrons /Xe+ = ~Me- = is the 

temperature of the electron gas. 

It is worth to notice, that by expressing factor ^ by 
the neutrino (antineutrino) energy rather than electron 
(positron) energy, we have just one formula, since both 
signs of AQ -I- me are covered, and f^, > 0. 

The neutrino spectrum from decay is proportional 
to: 



£l{^Q-£,)^{£,~ AQf -rr 
1 -f exp {£^ -AQ + ii)/kT 



-e{AQ-me-£^) 
(5) 

Figure[l]compares neutrino spectrum given by formula 
(|4]) with the more elaborated result of ^4^] for solar neu- 
trinos. Results are in good qualitative agreement. In 
both cases shown in Fig. [T] the neutrino spectrum is sim- 
ply a line of negligible ( Fig. [21 upper-left ) width. The 
horizontal axis in Fig. [T] is the difference between the Q- 
value (including me) and the neutrino energy, in keV. 
This is because for solar conditions the Q-value for ^Be 
capture is by many orders of magnitude larger than the 
temperature and the chemical potential of the electron 
gas. If we put ^Be into a plasma where kT or /i is com- 
parable to the Q-value, both capture rate and neutrino 
spectrum change dramatically, cf. Fig [21 In general. 




£„-AQ-m, [keV] 



FIG. 1: The normalized neutrino spectrum for solar ^Be elec- 
tron capture neutrino line, computed according to ^ (solid 
line) and the state-of-art result computed by Bahcall ( [48| . 
Eq. (46) ) is shown by dashed line. 



spectrum shape is a result of the competition between the 
Fermi-Dirac distribution and the unit step function Q in 
([3]). While e~ kinetic energy always adds to the neutrino 
energy, for low temperatures it is negligible compared 
to AQ + me- If the temperature becomes non-negligible 
compared to Q-values, say kT > 0.1 MeV, the thermal 
broadening due to kinetic energy of electrons becomes 
important and the capture rate is enhanced, cf. Fig. [21 
upper-right panel. For some of the laboratory stable nu- 
clei the electron (positron) capture might be possible for 
high energy electrons (positrons), from thermal distribu- 
tion tail. 

The increase of density result in large /ie, and that 
leads to a more visible effect. This is because most of the 
electrons, not just a small fraction from the tail, have 
large energies. The neutrino spectrum (Fig. [21 lower- 
left) has a very characteristic shape in this case, with 
sharp edge on the high £1, end. With the increasing /te 
progressively more nuclei become unstable to the elec- 
tron capture with the continuously growing capture rate. 
Lower-right panel in Fig. [H shows effect of large kT and 
/i. 

Anyway, possibly the most striking feature of Fig. [21 is 
not the shape of the spectrum but the dramating scale 
change on the vertical axis. Weak rates are extremely 
sensitive to both kT and /i, mainly due to phase-space 
factors dlini). 

In order to get combined NSE spectrum we have to 
sum up all terms (|4I5P for all relevant pairs of excited 
states, multiply them by the partition function and ma- 
trix elements, and then substitute into Eq. m with Xk 
obtained from NSE [42] . A typical behavior of the NSE 
i/e and De emissivities|38l l39l. l40l liol . [sol [Slj as a func- 
tion of Ye is presented in Fig. [3l As Ye decrease, electron 
neutrino flux (produced mainly in electron captures on 
protons and by heavy nuclei) also tends to decrease. On 
the other hand, a decrease in Ye cause an increase in 
the flux of S'e's. Usually antineutrino emissivity peaks 
due to beta decays of heavy nuclei and rise due to the 
positron capture on neutrons and by neutron decay, cf. 
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FIG. 2: The influence of tfie degeneracy (large n) and high temperatures (large kT) on the electron capture neutrino spectrum. 
Upper-left figure is for solar neutrinos (in laboratory conditions) and lower-left for cold degenerate electron gas. The upper-right 
figure refers to high temperature and the lower-right figure includes both degeneracy and high temperature. 



kT=0.6MeV, lgp=10 [g/cm^] 




0.0 0.2 0.4 0.6 0.8 

Ye 

FIG. 3: "(Color online)" Neutrino and antineutrino emis- 
sivities as functions of Y^. Critical Y^, defined as Ye — 0, is 
seen at the crossing point of the neutrino and the antineutrino 
particle emission rates. 



the neutrino astronomy, because they might lead to the 
strong neutrino and antineutrino emission without fur- 
ther neutronization, in agreement with constraints from 
nucleosynthesis. The increase of the flux (with de- 
creasing Ye) stops neutronization a little bit earlier than 
derived from e.g. the expansion of matter and the re- 
lated decrease in rates alone. The neutronization can 
also stop if Ye becomes too low and positron captures//?" 
decays start to dominate. Surprisingly, these critical Ye 
values (defined as Ye for which Ay_. = Ap^ , Fig. [3] and 
Eq. ^) vary in a broad range (Fig. 2]), reaching values 
close to Ye = 0.875 (primoidal BBN mixture of hydrogen 
and helium) for low densities and kT = 0.5 .. .0.8. On 
the other hand, for highest densities ( p > 10^^ g/cm^) 
and temperatures kT~0.8 MeV an equilibrium sets at 
Ye — 0.2 . . . 0.3. It is important to notice, that due to the 
low accuracy of the weak rates derived from FFN tables 
and the variability of the NSE state with Ye , Figure [4] 
provides only a very approximate outlook of critical val- 
ues. The critical value* is also very important for NSE 



Fig. [31 For almost all pairs {kT, p) we can find the value 
of Ye (Fig. m where the flux of Ve is equal to the flux of 
Ve- These threshold values are particularly interesting for 



* The state with Ye = is frequently refereed to as kinetic beta 
equilibrium. 
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TABLE I: Examples of neutrino and antineutrino spectra 



0.5 0.6 0.7 
kT [MeV] 



FIG. 4: "(Color online)" Critical (see Fig. |3] for explana- 
tion) for a range of considered temperatures and densities. 



timescales, as "stalled" Ye provide additional time with- 
out breaking assumption on the quasistatic Yg evolution. 

The competition between and v^. emission (Fig. 2]) 
is usually described in terms of the balance between 
electron captures (mainly on protons) and I3~ decays of 
the heavy nuclei [5^. However, for Ye outside range of 
0.35.. 0.45 the most important process leading to the Ve 
emission is the positron capture by neutrons. 



III. SPECTRA UNDER ASTROPHYSICAL 
CONDITIONS OF INTEREST 

We are able to compute approximate neu- 
trino/ antineutrino spectra for a wide range of as- 
trophysical phenomena if the NSE timescale is short 
compared to dynamic and weak timescales. Main 
limitation of our method is the neutrino trapping. The 
core-collapse supernovae and related phenomena, e.g. 
long gamma-ray bursts, are examples of objects where 
neutrino trapping is essential. We can use our method 
for initial infall stage of the collapse only. But as long 
as we are inthe free streaming regime this is the method 
of choice. We can produce much more detailed and 
accurate neutrino spectra (than hydrodynamic simu- 
lations itself) via postprocessing. The latter is trivial 
to parallelize, and allow to achieve greater accuracy- 
Our method can be applied to cosmological-like 1| 
neutrinos (Fig. [5]), the center of the pre-supernova 4| 
star (Fig. [S]) and typical conditions during type la 
thermonuclear explosions (Fig. [71 [5]). Other examples. 
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"Note: this example pushes our 
bility. More realistic spectrum is 
trapped and they begin to diffuse 



method to the Umits of applica- 
different, because neutrinos are 
rather than escape freely. 



less interesting since there are many known results 
Ji 111 01, 53, 54 55^ 56, 57, 58, 59, 60, 61, 62, 63, 64, 

mmmmMMM,MMMMM [Z3, Izi, m\ 

are provided for the core-collapse SN (Fig. [11]) and 
the merger of NS (Fig. The spectrum might be 

calculated for exotic conditions which are not related 
to any recently considered model as well, cf. Fig. [TOl 
Examples are listed in Table [J 

Illustrative example is provided using cosmological 
weak freezout values as input. Following we put 
Tg ~ 9.9, p = 0.008 g/cm^ and = 0.82. Neutrino 
and antineutrino spectrum in Fig. [5| is produced mainly 
from pair annihilation process. Therefore, spectrum is 
almost purely thermal. This is not a surprise, because 
thermal spectrum is what is expected for Big Bang neu- 
trinos. Thus, our method is working qualitatively well 
even in this extremal example. 

Pre-supernova stars are neutrino sources of particular 
interest [83|. The neutrino spectrum has been obtained 
using values at the center of a star during maximum com- 
pression stage. This stage is achieved just prior to the 
shell Si ignition above iron core, few hours before the 
start of the collapse. Spectrum for kT = 0.43 MeV, 
p ~ 7 X 10^ g/cm^ and Ye = 0.445 is presented in Fig. [Sj 
The important reference for these numbers is [3| ■ A strik- 
ing feature in Fig. [Sj is the significant contribution from 
heavy nuclei for both Ve and Pe spectra. This is partic- 
ularly important for the detection of these neutrinos, as 
previous studies [11] were based solely on the thermal 
emission. For the pair annihilation process dominates 
the high energy tail (fp^ > 10 MeV), and the number 
of detectable inverse-/? events in a standard large water 
Cherenkov detectors [H, [H, [sH] does not change. For 
a different detector design, e.g. a liquid scintillator [s^ 
the threshold might be as low as 0.2 MeV jSTj and the 
number of events will be much larger than anticipated 
from thermal processes only. The situation is even more 
pronounced for electron neutrinos. A large number of 
nuclei participate in massive electron captures leading 
to the flux that is two orders of magnitude larger than 
from pair process even for > 5 MeV. Therefore, the 
detection of fe's, previously rejected from analysis due 
to experimental difflculties, should be reconsidered. The 
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FIG. 5: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for BigBang nucleosynthesis era. 



contribution from free nucleons can be neglected in pre- 
supernova case. 

Another very important example of application for our 
method is the type la thermonuclear supernova. Two 
important regimes for thermonuclear burning in type la 
supernovae are deflagration and detonation. For defla- 
gration we use kT = 0.53 MeV and p = 2 x 10^ g/cm"^. 
For detonation a similar temperature of kT — 0.52 MeV 
has been used, but the density has been reduced due to 
pre-expansion to p = 7.8 x 10^ g/cm^. Ye — 0.5 was used, 
but it is important to point out that small neutroniza- 
tion is inevitable, because the ly^ flux dominates over the 
flux of De, cf. Fig. [8l Neutrino emission is very sensitive 
to these small changes. Model yl2 and n7dlrl0tl5c of 
[sol . [sH are the sources of the values used above. For 
a type la supernovae free nucleons are among the top 
neutrino sources, see Figs. [3 [H Two presented cases are 
related to detonation and deflagration. A lower density 
used for detonation stage (Fig. [7|) is the result of the 
white dwarf pre-expansion due to the previous deflagra- 
tion stage [80|- The high density in Fig.[8]is connected to 
the initial stage of subsonic nuclear burning in the pure 
deflagration model of [sOl. Only three nuclei contribute 
signiflcantly to the i^e spectrum in both cases: ^^Co, ^^Ni 
and protons. Relative contributions and total flux are dif- 
ferent, however. The neutrino flux per unit volume is four 
orders of magnitude larger for deflagration compared to 
detonation. Nevertheless, the deflagration involve tiny 
volume of the white dwarf only, while the detonation 
wave usually traverses entire star. Integrated flux might 
be similar, but this is model-dependent. Antineutrino 
spectrum is dominated by thermal processes: pair an- 
nihilation during detonation and plasnion decay during 
deflagration. The total flux is much smaller than for neu- 
trinos, and this imbalance causes to decrease. There- 
fore results from Figs. [7] and [H with assumed Ye = 0.5, 
should be taken with care. For example, the NSE abun- 
dance of ^^Co drops rapidly in the range Yg = 0.5 . . . 0.47. 
A more detailed investigation of type la neutrinos shows 



also an important contribution from free neutrons to the 
anti-neutrino spectrum above 2 MeV. 

Now, we study neutrino spectrum for the accretion disk 
formed in NS-NS merger. Needed data are taken from 
Fig. 1 of [Hi: temperature of kT = 1.0 MeV, p = 10^° 
g/cm'^ and Ye = 0.05. Similar results are expected for the 
neutron star - black hole mergers and other phenomena 
forming low Yg, dense, high temperature accretion disks. 
The neutrino spectrum is a result of pair process and 
electron captures on protons. The antineutrino spectrum 
is heavily dominated by the neutron decay and positron 
captures on neutrons. The gap 0.8 < < 1.8 MeV 
is filled by processes involving heavy nuclei. Moreover, 
antineutrino flux is much larger compared to the neutrino 
flux. The spectrum peaks at ~ 5 MeV, providing 
interesting candidate for the neutrino detection using the 
inverse /? decay. 

Another example (Fig. fTU]) . not related to a particular 
astrophysical phenomena, shows the importance of ther- 
mal processes and those involving free nucleons. The 
antineutrino spectrum, especially the high energy end 
due to the positron capture is particularly important. 
Spectral features of this process should be interesting for 
future neutrino astronomy, based on gigantic De water- 
based detectors |88i, HE (IS ■ 

The core-collapse process is poorly described using our 
method, but we have provided an example for the sake 
of the completeness. The calculated neutrino spectrum 
in Fig. [11] has a complex multi-peak structure. This is 
in contrast to the results of the more sophisticated neu- 
trino radiation transport results, which are always single- 
peaked. This can be explained by: (1) smoothing na- 
ture of the diffusive transport and (2) too small energy 
resolution (to few energy bins) of the transport codes 
used in simulations. The high energy neutrinos seen in 
Fig. [Tl] are in reality downscattered to much smaller en- 
ergies. The same applies to antineutrinos. Additionally, 
ly-D pairs are created in the process of collisions between 
neutrinos and electrons, and between pairs of neutrinos 
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FIG. 6: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for presupernova star. 




FIG. 7: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for detonation stage of thermonuclear supernova in delayed-detonation class of models. 



[9ll |. This leads to the energy exchange between flavors, 
and realistic Vx spectra^ are not as distinct as those from 
Figs.[TlJ Factors that block outgoing neutrinos and could 
shape the neutrino spectrum under such extreme condi- 
tions were omitted. Clearly, our method is not working 
for the core-collapse supernovae, as anticipated. 



IV. CONCLUSIONS 

One of our important conclusions is related to typi- 
cal way of publishing data on weak nuclear processes in 
astrophysics. This approach dates back into year 1980, 
and was introduced in the famous paper [3a]. Tables 
published by the FFN become standard in modern astro- 
physics. Upgrades [SS*, "4^, "oS] did not change structure of 
FFN tables. Unfortunately, FFN grid using mere 13x11 



^ Muon and tau spectra are almost identical to the thermal elec- 
tron flavor spectra, except for smaller integrated flux. 



points is not enough to obtain precise results, as noted al- 
ready by the FFN authors [4l| . While we understand rea- 
sons to preserve this standard for 30 years, "reverse en- 
gineering" of FFN-like tables to get spectrum, as well as 
complicated interpolating procedure is impractical now. 
If one wants to calculate the spectrum precisely, with- 
out analytical approximate formula for individual nuclei, 
pre-calculated tables are useless. Much more convenient 
is the following set of data: 

1. energy and spins for ground and excited states 

2. weak transition matrix elements between all rele- 
vant pairs of the excited states for the parent and 
daughter nuclei 

Alternatively, tabulated spectrum for all T — pY,, pairs 
would be a good choice, with amount of stored data up 
to several megabytes. While such approach will increase 
amount of published numerical data by a factor of ~ 10, 
it would remove any ambiguity in the representation of 
the spectra. 
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FIG. 8: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for early deflagration stage of thermonuclear supernova. 
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FIG. 9: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for neutron star merger. 



The inspection of virtually any of the figures presented 
here (Figs.l MTT]) clearly show the importance of both nu- 
clear and thermal processes. The thermal emission and 
captures on free nucleons and nuclei should be included 
in consistent calculations. However, depending on the 
subject, all combinations of these can be found in astro- 
physical applications. For example, type la supernova 
simulations include NSE emission but older simulations 
neglect neutrino emission at all or include electron cap- 
tures only [23| . Other important regimes, core-collapse 
and pre-supernovae frequently neglect positron captures, 
particularly on neutrons. Estimates of the neutrino sig- 
nal in detectors from pre-supernovae rely purely on ther- 
mal emission 83,, ,83, . 

The ultimate goal which is beyond scope of the arti- 
cle is to know exactly (not approximately!) the neutrino 
spectrum from weak nuclear processes under NSE. In the 
past weak rates were usually integrated and only the total 
neutrino flux (particles and energy) has been tabulated 
and presented to the public. We argue again, that this 
is not the best approach if one wants to calculate the 
neutrino spectrum. Without full input used to calculate 



weak rates we are unable to restore information lost in 
the integration. Typical (FFN-like) weak interaction ta- 
bles are not sufficient. Tables of the excited states, spins 
and weak matrix elements for all considered nuclei will 
allow researchers to calculate both neutrino/ antineutrino 
spectra and customized weak interaction rate tables. 

Weak rates prepared in the FFN fashion (i.e. all pub- 
lished rates [H, [3^, [13]), even those with tabulated ef- 
fective Ig (ft) , do not facilitate estimates of the neutrino 
spectrum. This is not surprising, because these rates 
were prepared for a different purpose: the neutrino en- 
ergy loss and neutronization. Maximal information on 
the spectrum extracted from FFN-like tables can be ex- 
tracted as described in the paper accompanying paper 
[43 |. We re-tabulate effective Ig (ft) values and effective 
Qef f-yalues for every grid point to get from (|3]) or ([5]) the 
original total rate and average neutrino energy. If the to- 
tal rate is not dominated by the captures we switch from 
dH) to ([5]) . This approach produces significant side effects 
if capture and decay rates are comparable. The neutron 
provides good example. Due to the non-negligible con- 
tribution of Pe's from the neutron decay, the average en- 
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FIG. 10: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions of large 
Ye and high temperature. 




FIG. 11: "(Color online)" Neutrino (left) and antineutrino (right) spectrum emitted per unit volume under conditions typical 
for the infall stage of the core-collapse. 



ergy differ from that deduced from pure positron capture. 
Therefore the effective spectrum has a variable effective 
Q-value. The realistic positron (and electron as well) 
capture spectrum always starts with energy equal to the 
lowest Q-value. To sum up, the obvious next step in the 
research is to give up pre-calculated tables of weak rates 
and to re-calculate the neutrino spectrum from scratch, 
using nuclear data and weak matrix elements as an input. 

Despite these difficulties, we obtained new results. 

1. we get interpolating procedures for NSE abun- 
dances with number of convenient features: the 
ability to pick out of NSE selected nuclei, the 
computational time scaling linearly with the num- 
ber of nuclides and independent of the position in 
T - p-Ye space for full Y^ = 0.05 . . . 0.95 range 




2. the energy spectrum, fluxes, mean energies etc. 
of the emitted neutrinos and antineutrinos sepa- 
rately for i/e and I'e 



Our analysis was meant to be general, but we can 
identify some possible astrophysical targets for presented 
methods. The NSE neutrino spectrum would be a good 
approximation for massive stars after Si burning and 
thermonuclear supernovae. A related research is under- 
way. Procedures developed here will be useful for the 
analysis of neutrino signals from X-ray flashes, neutron 
stars, merger events, accretion disks and some types of 
cosmic explosions, e.g. pair-instability supernovae. 

The electron antineutrino emission due to the positron 
capture on neutrons provides strong and relatively high- 
energy flux for surprisingly large volume in kT — p — Ye 
space. Needed thermodynamic conditions: kT > 0.6 
and p > 10^ g/cm^ can be met in many astrophysi- 
cal objects. Megaton-scale neutrino detectors [93] will 
search for antineutrinos with energy > 1.8 MeV. 
The detection of strong I'e-flux above 5 MeV produced 
mainly by captures on protons and heavy nuclei is stan- 
dard in water Cherenkov [H, [s^, [sO] or liquid scintillator 
[H, [H, detectors. Therefore further investigation of 
NSE neutrinos, particularly in the unexplored region of 
large 0.87 > Ye 3> 0.55 should give researchers some 
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additional hints for the existence (or non-existence) of 
detectable astrophysical antineutrino sources. 
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